Non-linear photonic crystals as a source of entangled photons 
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Non-linear photonic crystals can be used to provide phase-matching for frequency conversion in 
optically isotropic materials. The phase-matching mechanism proposed here is a combination of form 
birefringence and phase velocity dispersion in a periodic structure. Since the phase-matching relies 
on the geometry of the photonic crystal, it becomes possible to use highly non-linear materials. This 
is illustrated considering a one-dimensional periodic Alo.4Gao.6As / air structure for the generation 
of 1.5 fim light. We show that phase-matching conditions used in schemes to create entangled 
photon pairs can be achieved in photonic crystals. 
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Polarization entangled photon pairs play a central role 
in testing the foundations of quantum mechanics and in 
implementing quantum information protocols 0. They 
are likely to remain an appealing resource for practical 
quantum information science since they interact very lit- 
tle with the environment, propagate easily over long dis- 
tances and since photon polarization is easily manipu- 
lated in experiments. A popular method to create such 
pairs is by down-conversion in non-linear crystals. For 
this process to be efficient the down-converted photons 
generated by pump photons at different points in the 
crystal have to be in phase with each other. This leads 
to the "phase-matching" condition kp — ki + k2, for the 
wavevectors of the pump (kp) and downconverted pho- 
tons (fci , ^2) . In general this condition will not hold owing 
to normal index dispersion which makes \kp\ > |fci|-l-|fc2|- 
The current schemes for down-conversion employ the nat- 
ural birefringence of specific non-linear crystals, like /3- 
Barium-Borate (BBO), to compensate for this effect. 

In one of these schemes 0, Q , down-converted photons 
emerging in a particular pair of directions are entangled 
in polarization: 

\^) = ^i\V)^\H)2+c^'>\H),\Vh), (1) 

where V and H denote the polarization state of particles 
1 and 2 and is a phase factor that can be set by passing 
one of the photons through a phase-plate. 

Although birefringent non-linear crystals have proven 
to be a successful source in many proof of principle ex- 
periments the demands of both practical quantum infor- 
mation and schemes for the implementation of all op- 
tical quantum computing 0| make it desirable to pro- 
duce sources that are both more efficient and more eas- 
ily integrable on, for example, an optical chip. The re- 
quirement that both the birefringence and the x^^^ non- 
linearity be naturally present severely limits the possibil- 
ities for improvement of sources based on existing crys- 
tals. Semiconductors such as GaAs or GaP have a 
typically 200 pm/V IBI, about two orders of magnitude 



larger than that of commonly used crystals such as BBO, 
2.2 pm/V 0. This, together with the existence of well 
developed micro-fabrication techniques for these materi- 
als, makes it attractive to explore ways of creating semi- 
conductor based entangled photon sources. In these ma- 
terials, which have no natural birefringence, the condi- 
tions for phase-matching must be created artificially. In 
this letter we show how photonic crystals can be used to 
achieve this. 

Photonic crystals are materials with a periodic varia- 
tion in refractive index on the scale of the optical wave- 
length. In layered structures the different boundary con- 
ditions at the interfaces for the two polarizations leads 
to an effective or "form" birefringence. The application 
of this effect to phase-matching has been explored in the 
limit that the optical wavelength is much larger than the 
periodicity 0, H, ■ When the wavelength is compara- 
ble to the periodicity it was suggested that a reciprocal 
lattice vector can be added to the phase-matching con- 
ditions Alternatively a similar relation was found 
for periodically poled structures that have a periodicity 
in x'"^\ leading to quasi-phase-matching conditions 
Our work proposes how a combination of form birefrin- 
gence and the strongly altered phase velocity in a pho- 
tonic crystal can be used to phase-match down-converted 
light and generate entangled photons. 

We focus on one-dimensional structures for sake of sim- 
plicity, although the concepts can be extended to two- 
and three-dimensional photonic crystals. The periodic 
variation in refractive index leads to Bragg scattering of 
the light and wave propagation becomes best described in 
terms of a photonic bandstructure. If the wavelength is 
comparable to the periodicity of the structure, the prop- 
agation of light is strongly affected, leading to the exis- 
tence of a range of frequencies, known as a stopband, for 
which light does not propagate. 

The propagation of electromagnetic waves in a pho- 
tonic crystal is described by Maxwell's equations with a 
periodic dielectric function e(r). The general solution is 
given in terms of Bloch waves labeled by a frequency uj 
and a Bloch-wave vector K and consists of a plane wave 



multiplied by a function that has the periodicity of the 
photonic lattice. The Bloch wave can be written as a 
Fourier sum over the reciprocal lattice vectors G: 

Egir, t) = e-(^-^-*) J2 ^K^G^-''^ ^ (2) 

G 

where e^^ g are the Fourier coefficients for each of the 

space harmonics. For a one-dimensional structure K and 
the coefficients g can be obtained by finding the eigen- 
values of a 2 X 2 transfer matrix for each polarization 
separately 0. The effect of index dispersion, essential in 
any discussion of phase-matching, can be easily incorpo- 
rated as the transfer matrix is defined at each frequency 
separately. The obtained wavevector can be decomposed 
in components Kz perpendicular and k\\ parallel to the 
layers of the photonic crystal. The dispersion relation for 
transverse electric (TE) or ordinary (o) polarized waves, 
that have the electrical field in the plane of the layers, is 
given by 

Kz{ijJ,k\\) — — arccos[cos(fci2a + /c22&) ~ 
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where kiz and k2z are the plane wavevector components 
perpendicular to the interfaces in medium one (of refrac- 
tive index rii) and medium two (of refractive index 712): 

An analogous expression can be derived for the orthog- 
onal transverse magnetic (TM) or extraordinary (e) po- 
larization, that results in different propagation constants. 
This important polarization dependence arises from the 
different boundary conditions at the interfaces for the two 
polarizations and will play a key role in the generation 
of polarization entangled photon pairs. 

Figure^shows a bandstructure derived from Eq.Elfor a 
periodic structure with alternating layers of Alo.4Gao.6As 
and air. The fill fraction of Alo.4Gao.6As is 0.656. The 
gray area corresponds to propagating solutions for TM 
waves (left panel, negative fc||) and TE waves (right panel, 
positive fc||). They coincide for normal incidence (fcy — 0). 
The solid lines in the figure correspond light in vacuum 
[lo = c\k\) and divide the modes in those accessible to 
waves from outside the crystal, and those that are con- 
fined by total internal reflection. The latter modes can 
be accessed from the side of a sample or by using a set 
of prisms [l^ . 

The frequency is specified in multiples of ttc/A, where 
A is the periodicity of the structure. The structures were 
designed by first considering the refractive index to be 
equal to that at the pump frequency ojp {n = 3.4 at 
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FIG. 1: Bandstructure of a periodic structure with 
Alo.4Gao.6As and air layers for TM (left) and TE (right) po- 
larized light. The grey area corresponds to propagating so- 
lutions at a given frequency uj and wavevector fcy . The solid 
lines indicate the "light line" in vacuum: uj = c\k\. 

A = 750 nm). This allows the selection of a period- 
icity A such that the pump photons propagate in the 
structure. Once the periodicity is fixed, the frequency 
dependent refractive index is used to calculate disper- 
sion surfaces to search for phase-matching conditions at 
relevant frequencies. The choice of material is motivated 
by the large second order non-linearity together with the 
fact that it is transparent at a pump wavelength of 750 
nm, to allow degenerate down-conversion in the impor- 
tant telecommunication window around 1500 nm. More- 
over, the fabrication of such crystals seems to be feasible 
by starting with a periodic structure that has layers of 
AlAs and Al o 4G ao.6As followed by selective wet etching 
of the AlAs 13 . The fiU fraction was chosen in order 
to maximize birefringence in the long wavelength limit 
(A ^ A) , where the Bloch waves are plane waves and ex- 
perience an effective medium that behaves as a uniaxial 
birefringent material. 

The leading terms in the Bloch wave expansion in Eq.[21 
will phase-match when Kp + Gp = Ki + Gi + K2 + G2 
where the Gs correspond to the leading eg ^s, ensuring 
efficient down-conversion into these modes. Two regimes 
are explored that lead to the generation of entangled 
photon pairs. We first discuss the long wavelength limit 
(A ^ A) followed by the case where A ~ A. 

Figure [3 shows dispersion surfaces in the long wave- 
length limit for a structure with periodicity A = 18.75 nm 
and a pump wavelength of 750 nm (0.05 ttc/A). One 
set of dispersion surfaces for the down-converted light, 
at 1500 nm wavelength, is drawn at the origin, another 
is drawn at the end of the pump wavevector in order 
to obtain a simple geometric construction of the phase- 
matching condition. Wavevectors satisfying the phase- 
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FIG. 2: Dispersion surfaces for degenerate down-conversion 
at 1500 nm in a one-dimensional Alo.4Gao.6As / air photonic 
crystal with a periodicity A = 18.75 nm and a fill fraction of 
AlQ.4Gao.6As of 65.6%. Dispersion surfaces are shown for the 
pump at a wavelength of 750 nm (thick lines) and for degener- 
ate down-conversion at 1500 nm (thin lines) . The dashed lines 
correspond to TE (ordinary) polarization, while the solid lines 
correspond to TM (extraordinary) polarization. The pump 
wavevector indicated by the arrow is extraordinary. 



matching condition correspond, in such a diagram, to 
the intersections between the dispersion surfaces of the 
down-converted photons. 

The intersections of the thin dotted hues represent 
the wavevectors of two ordinary (o) polarized photons, 
created by the extraordinary (e) polarized pump field. 
This situation is usually referred to as non-colinear type-I 
down-conversion. Since the structure is symmetric under 
rotations about the optical axis parallel to K^, the cir- 
cles and ellipses in Fig. [3 should be seen as spheres and 
ellipsoids. The intersections of the two down-conversion 
spheres for o polarization will give a circle, represent- 
ing a cone of down-converted light centered on the pump 
wavevector. Every pair of photons is emitted in such a 
way that they obey the phase-matching condition and so 
are always diametrically opposite each other about the 
pump. 

For noncolinear type-II down-conversion, the pump 
field is e polarized and the down-converted photon pairs 
each have one o polarized and one e polarized photon. 
This case corresponds to the intersections of the dashed 
thin lines with the solid thin lines in Fig. [3 There are two 
pairs of such intersections, each representing an intersec- 
tion between a sphere and an ellipsoid. They define two 
slightly distorted cones of down-converted light whose 
central axes are not aligned. For the direct generation of 
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FIG. 3: Dispersion diagram in the repeated zone scheme for 
a pump wavelength of 750 nm. The periodicity of the crys- 
tal is A = 187.5 nm and the fill fraction of AIo.4Gao.6As is 
a = 0.656. The pump wavevector is indicated by the arrow. 



polarization entangled photons as in Eq. ^ the e and o 
photons in a pair must be emitted in directions in which 
e photons and o photons are indistinguishable except by 
their polarization. This happens at the intersections of 
the two distorted cones discussed above. 

Figure 13 shows dispersion surfaces for a photonic crys- 
tal with a periodicity A = 187.5 nm, a factor of 10 greater 
than that of Fig.|21and comparable to the wavelength of 
the pump light. The frequency of the pump w = 0.5 
TTc/A (A = 750 nm) is now above the first stopband for 
normal incidence. The dispersion surfaces for the pump 
are no longer continuous as some directions are excluded 
by Bragg reflection. 

For the downconverted photons at a wavelength of 
1500 nm, the Bloch waves are essentially plane waves 
with wave vectors Ki and K2 ■ The phase-matching con- 
dition in this case reduces to: Kp + Gp = Ki + K2- The 
pump wavevector, indicated by the arrow is chosen to 
end in the second Brillouin zone in order to select the re- 
ciprocal lattice vector that corresponds to the dominant 
space harmonic of the Bloch wave in Eq.[5] In principle, 
phase-matching can also be achieved by choosing another 
reciprocal lattice vector, but the amplitude of the down- 
converted light will be reduced. 

The dispersion surfaces of down-converted photons in 
Fig. 121 indicated by the thin lines resemble those of the 
long wavelength limit. Again there are two different 
types of intersections, corresponding to Type-I (e ^ 00) 
and Type-II (e eo) down-conversion. As discussed 
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in the long wavelength limit, polarization entangled pho- 
tons may be collected at the intersections of the distorted 
cones. 

In the long wavelength limit the phase-matching relies 
solely on form birefringence. However, when the peri- 
odicity of the crystal is comparable to the wavelength of 
the pump, the phase-matching also relies on a substantial 
change in phase velocity. For frequencies below the stop- 
band the phase velocity is lower than that in the long 
wavelength limit, which works against phase-matching 
since it adds to the effect of normal index dispersion. 
However, if the frequency is above the frequency of the 
stopband the phase velocity is increased and compensates 
the effect of index dispersion. For large enough index con- 
trast this can be used to achieve phase-matching where 
all photons have the same polarization, i.e. e — > ee or 
o oo. 

We will now estimate the efficiency of our scheme com- 
pared to that of the scheme based on BBO. There are 
four factors that contribute to the relative amplitude: (i) 
The ratio of the value of the x'^^-''s. (ii) The efficiency 
per unit length should be multiplied by the fill fraction 
of x^^-' material (0.65). (iii) The scheme phase matches 
the leading term of the Bloch waves. The magnitude of 
this term is obtained by Fourier analysis of the Bloch 
wave. This gives a factor of ~1 and a factor ~0.21 for 
the structures in Fig. Inland Fig. O respectively, (iv) As 
for any non-linear material we need to consider the tensor 
properties of the non-linearity. Alo.4Gao.6As has the 43to 
point group symmetry and for the conventional (100) sur- 
face orientation and our directions this gives a factor of 
'--^O.SS and ~0.53. Multiplying all factors and squaring 
to get the relative efficiency gives ^--^2500 and ^^50 times 
the efficiency for BBO. 

The example structures we used to illustrate our 
scheme have by no means been tailored to maximize the 
efficiency. The large value of x*'^^ in Alo.4Gao.6As leaves 
ample room for tailoring the geometry of the structure 
to suit a particular application, whilst keeping the pro- 
cess efficient. It should be noted that there is a trade- 
off between increasing the strength of the interaction of 
light with the photonic crystal and the loss in efficiency 
associated with the undesired terms in the Bloch wave 
expansion (c.f. the efhciency for Fig. [21 and Fig. |3J. An 
interesting extension of our work is to consider downcon- 
version of pulsed gaussian beams, which is accompanied 
by spatio-temporal walk-off effects that affect the qual- 
ity of entanglement and the efficiency of the process. It 
seems promising to explore the difference between group 
velocity and phase velocity in a photonic crystal to re- 
duce these effects. 

The concepts introduced in this letter can be extended 
to two- and three-dimensional photonic crystals. The 
change in phase velocity is a result of Bragg reflection 
common to all photonic crystals. Form birefringence re- 
lies on polarization dependent boundary conditions and 



is observed in one and two-dimensional photonic crys- 
tals. Form birefringence can be introduced in three- 
dimensional structures by making the individual scatter- 
ers elongated which can be achieved in microfabricated 
structures J4j or by deforming a colloidal crystals jlS] . 

We presented a phase-matching scheme that uses the 
combined effects of form birefringence and a change in 
phase velocity in photonic crystals to enable the produc- 
tion of entangled photon pairs. Both effects can be tuned 
by designing an appropriate structure and are decoupled 
from the intrinsic non-linearity that is determined by the 
choice of the constituent materials. We focused primar- 
ily on III-V semiconductors and on AlGaAs in particu- 
lar. These materials have a high non-linearity together 
with a large linear refractive index which contribute to 
a strong interaction with light. Microfabrication tech- 
niques are well developed for these materials and an ex- 
perimental realization seems feasible. If successful, the 
structure could be implemented on optical chips as an 
integrated source of entangled photons. 
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